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^ . Abstract 

We study the octet SS — PP correlator within resonance chiral theory up to the one- 
loop level, i.e., up to next-to-leading order in the 1/Nc expansion. We will require that 
our correlator follows the power behaviour prescribed by the operator product expansion at 
high euclidian momentum. Nevertheless, we will not make use of short-distance constraints 
from other observables. Likewise, the high-energy behaviour will be demanded for the whole 
correlator, not for individual absorptive channels. The amplitude is progressively improved 
by considering more and more complicated operators in the hadronic lagrangian. Matching 
the resonance chiral theory result with chiral perturbation theory at low energies produces 
the estimates L 8 (u) su ^ = (1.0±0.4) • 10~ 3 and C 38 (u) su ^ = (8±5)-10" 6 for u = 770 MeV. 
The effect of alternative renormalization schemes is also discussed in the article. 
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1 Introduction 

The effective field theory (EFT) approach is a very powerful tool for the investigation of Quantum 
Chromodynamics (QCD) at long distances. Chiral Perturbation theory (%PT) [U [2], [3] is the EFT 
for the description of the chiral (pseudo) Goldstones in the low energy domain p 2 <C A 2 H ~ 1 GeV 2 , 
with Ah typically the scale of the lowest resonance masses. The calculation of the QCD matrix 
elements is then organized at long distances in growing powers of the external momenta and light 
quark masses. Recent progress has allowed to carry up to 0(p 6 ), i.e., up to the two-loop 
level H El El Ej. 

In the intermediate resonance region, Ah < E < 2 GeV, ;\RT stops being valid and one must 
explicitly include the resonance fields in the Lagrangian description. Unfortunately, this is not a 
straightforward process because there is no natural expansion parameter in this region as several 
relevant mass scales appear in this range (resonance masses, momenta, widths, the characteristic 
%PT loop scale A x ~ AttF...). Resonance Chiral Theory (R%T) describes the interaction of reso- 
nance and pseudo- Goldstones within a general chiral invariant framework [HE]. Alternatively to 
the chiral counting, it uses the 1/Nc expansion of QCD in the limit of large number of colours [10] 
as a guideline to organize the perturbative expansion. At leading order (LO), just tree-level 
diagrams contribute while loop diagrams yield higher order effects. Integrating out the heavy res- 
onance states leaves at low energies the corresponding chiral invariant effective theory, %PT. Many 
works have investigated various aspects of R%T: equivalence of formalisms [91 [HI [121 EE] ; Green 
functions pH UM UM EH CES1 Ell [20] ; applications to phenomenology [III [2H [22l [23l [2H [25l [261 [27] ; 
determination of chiral low-energy constants (LECs) at NLO in 1/N C [2H EH EH ED E2]; de- 
termination of the one-loop ultraviolet divergence structures in the generating functional [S3] ; 
implications about the renormalizability [SUES]; possible issues with extra degrees of freedom in 
the renormalized propagator [36l [37] ; renormalization group studies [38J . 

The infinite tower of mesons contained in large- Nq QCD is often truncated to the lowest states 
in each channel, usually named as single resonance approximation (SRA). This approximation has 
led to successful predictions of 0{p A ) and 0{p & ) low-energy constants (LECs) [SI El Ell ESI EH]- 
However, the study of Regge models with an infinite number of mesons has shown that if one keeps 
just the lightest states with exactly the same couplings and masses of the full model then one finds 
problems in the short-distance matching and wrong values are obtained for the LECs [10] . Thus, 
in a high-energy matching with the operator product expansion (OPE) [H] the parameters of the 
truncated theory will be shifted in order to accommodate the right short-distance dependence. 
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Chiral symmetry ensures the proper low-momentum structure of the R%T amplitudes around 
p 2 = but their high energy behaviour is not fixed by symmetry alone. In that sense, the matched 
amplitude can be understood with the help of Pade approximants as an rational interpolator 
between the deep Euclidean p 2 = — oo and p 2 = [331 S3]- The Weinberg sum-rules (WSR) [32] 
yield the most convenient parameters for the interpolation rather than the accurate determinations 
of the resonance couplings. Furthermore, the R%T couplings for the lightest mesons are expected 
to be in better agreement, whereas the parameters from the highest excitations may lie far from 
their right values [33] ■ 

The connection of the R%T amplitudes with the operator product expansion (OPE) at high 
energies seems a priori a useful procedure to include extra information from QCD in the resonance 
theory. It allows to fix combinations of couplings (e.g., through WSR), decreasing the number of 
unknown parameters in the analysis. However, large- Nc QCD has an infinite number of hadrons 
and in order to reproduce the full large-iV^ theory one must consider the tree-level exchanges of 
heavier and heavier resonances. In the hadronical ansatz approach, one adds more and more poles 
to the rational approximant [431 [33] • Equivalently, this can be realized within the quantum field 
theory framework as a generating functional with a lagrangian including interaction operators 
J — Rj that couple the external current source J and heavier and heavier resonances Rj (e.g. of 
the form c m j( SjX+ ) for the SS correlator). 

The extension of R;\T beyond the tree level approximation still needs to be worked out in detail. 
Although some theoretical issues on the renormalizability of R;\T still need further clarification [331 
1351 14"5] . several chiral LECs have been already computed up to NLO in l/Nc through quantum 
field theory (QFT) one-loop calculations [291 [30] an d dispersion relations [311 [32]. In this article 
we will focus our attention on the chiral octet SS — PP correlator (for instance, with 1=1), 
which in the chiral limit is determined at low energies by the 0(p 4 ) and 0(p 6 ) LECs, respectively, 
by L s [3] and C 38 [3]. The correlator is computed up to next-to-leading order in l/Nc (NLO) and 
the chiral limit will be assumed all along the article. 

At the one- loop level -NLO in l/Nc-, one needs also to devise a procedure to reach the infi- 
nite resonance limit of large-iVc QCD. In the case of two-point Green-functions, the imaginary 
part of the one-loop diagrams is given through the optical theorem by the square of two-meson 
form-factors computed at tree-level. Thus, based on a dispersive approach, one may add the con- 
tribution to the spectral function from higher and higher two-meson absorptive cuts by providing 
the corresponding form-factors [311 [32] . This would be, in some sense, the natural extension of 
the minimal hadronical ansatz |44j to the one-loop situation. In a previous computation of the 
octet SS — PP correlator up to NLO in l/Nc, the intermediate two-meson channels were analyzed 
individually [3T]. The corresponding tree- level form-factors were made to vanish appropriately at 
high energies [321 145] . This allowed to recover the correlator from its spectral function through 
an unsubtracted dispersion relation. However, in general, it is not always possible to fulfill the 
high-energy constraints for all the form-factors at once 0- Only the two-meson absorptive cuts 
with at most one resonance (tttt and Rtt) were considered in Ref. [31], as the RR' channels have 
their thresholds at (M R + M R i) ~ 2 GeV and are suppressed at low energies. Likewise, the 

1 In the case of the scalar and pseudo-scalar form- factors, it is still possible to impose the right high-energy 
behaviour to all the form- factors if one considers operators with two and three resonance fields £rr> and £rrir<< [321 
136] . Nonetheless, there is no consistent set of constraints for all the vector and axial- vector form-factors if only 
a finite number of resonances is considered [351 H5] . A similar kind of inconsistences was found in the study of 
three-point Green- functions at large Nc [19j . 
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short-distance constraints from VV — AA Weinberg sum-rules and the 7T7T vector and the scalar 
form-factor were used there in order to fix some of the couplings appearing in the analysis. 

In the quantum field theory approach proposed in this work, one has a mesonic lagrangian 
which at the classical level generates the large-iVc amplitudes and whose quadratic fluctuations 
around the classical field configuration provide the one-loop corrections [33]. The complete QCD 
generating functional is approached as one adds more and more hadronic operators to the action. 
Eventually, one should add the infinite number of possible terms of the given 1/Nc order under 
consideration. For instance, the Sim interaction (provided by Cd{Su^u^) |8J) is of the same 
order as in 1/N C as the SPn vertex (given by the Af p ( {V^S 1 , P}u^ ) operator [TSJ E21 USD- 
Notice that one never has a complete description with a finite number of operators. The basic 
lagrangian Lq + £r with at most one resonance field in each term [8] provides an incomplete 
description of the Rw channels, as the possible diagrams with R' resonances exchanged in the s- 
channel are missing [311 [32] . This requires the incorporation of operators Crr> with two resonance 
fields [T51 E21 HE] • In the same way, the RR' absorptive cuts are now badly described without the 
Crrir" terms with three resonance fields. 

The chiral structure of the lagrangian ensures the right structure at long distances. On the 
other hand, we will impose that the correlator follows the short-distance behaviour prescribed 
by the OPE. The one-loop R/\T amplitude will be used as an improved interpolator between 
low and high energies. The resonance couplings become then interpolating parameters that must 
approach their actual values in the full QCD as more and more operators are added to the R/\T 
action. On the contrary to what was done in former works [311 [32], the short- distance matching 
will be carried out in the present article for the total correlator and spectral function [37], rather 
than for individual channels. Likewise, we will not use the short-distance constraints from other 
amplitudes to fix the couplings in the one-loop correlator. We will work within the SRA, including 
just the chiral Goldstones and the lightest multiplets of scalar, pseudo-scalar, vector and axial- 
vector resonances. In a first step, the SS — PP correlator will be computed at NLO in 1/Nc with 
the simplest R/\T lagrangian, with operators with at most one resonance field (Gy, c m , d m ...) [8]. 
This provides the proper structure for the intermediate tree- level exchanges (it, S, P one-particle 
channels) and the two-Goldstone cut nn. However, this simple lagrangian fails to describe the 
R-k and RR' channels as the lagrangian [S] makes their form-factors behave like a constant or 
like a growing power of the momentum at high energies [30l [3T| 1321 [3H1 EE]- This will be partly 
cured by the consideration of X RR operators with two resonance fields [T5l 1321 [3H1 [3H], which 
now allow an appropriate description of the Rw channels, though the RR' ones still behave badly. 
Although these cuts with two resonances were neglected in the dispersive approach [31] , removing 
part of the one-loop diagrams is not theoretically well defined and may lead to inconsistences in 
the renormalization of the QFT. Furthermore, it is not trivial that the effect of the RR' cuts in 
the short- distance matching is fully negligible. Hence, all the possible diagrams contributing to 
the correlator up to NLO will be kept in our study. 

The amplitude is first computed within the usual subtraction scheme of [2] (denoted 

for simplicity as MS all along the article). However, though equivalent at low energies, some 
appropriate schemes will be found more convenient: pole masses and other schemes that minimize 
the uncertainties derived from the short- distance constraints. This will help us to determine the 
0(p 4 ) and 0{p & ) LECs, respectively L 8 (fi) and C 38 (/u). The high-energy constraints and their 
meaning will be discussed and the convergence to full large-iV^ QCD will be tested as more and 
more hadronic operators are added to the R%T action. This work is thought as a complementary 
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and an alternative approach to the dispersive analysis in Ref. [3 lj . 

The article is organized as follows. Resonance chiral theory is introduced in detail in Sec. [2j 
The octet SS — PP correlator is defined in Sec. [3] and its one-loop R%T computation is provided 
in Sec. HI The high-energy constraints and low energy expansions are respectively given in Sees. 
and [6j The contributions from operators £rr> with two resonance fields have been singled out in 
Sec. [7] to ease the main argumentation of the article. Finally, the phenomenological analysis is 
given in Sec. [8] and the conclusions are provided in Sec. [9j Some technical results are relegated to 
the Appendices. 



2 Resonance chiral theory lagrangian 

Within the large-iVc approach the mesons will be classified within U(3) multiplets. The chiral 
Goldstone bosons are introduced by means of the basic building block, 
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with x = 2B (s + ip) containing the scalar and pseudo-scalar external sources, s and p respectively, 
the right and left sources r M and ^ providing the vector and axial-vector external sources, v M = 
i( r M_|_£M) anc i a M = I( r M_£A») respectively, and F^ U R the corresponding left and right field-strength 
tensors. 

The Goldstone bosons are parametrized by the elements u(<f>) of the coset space U(3)l x 
U(3)r/U(3)v, transforming as 



u{4>) h- VRU^Kg^Y 1 = hfafluifiV! 



R 



(4) 



under a general chiral rotation g = (Vl, Vr) C G in terms of the U{3)y compensator field h(g, • 
This makes the tensors X = x±, f± to transform covariantly in the form, 



X 



h(g,4>)Xh(g,4>r 1 . 



(5) 
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2.1 Leading order lagrangian 

For the classification of the vertices entering in the tree-level and one-loop amplitudes it will be 
useful to organize the operators of the R%T lagrangian according to the number of resonance 
fields: 

C = C G + C R + Crr> + .» (6) 

where Lq only contains Goldstone bosons and external sources, Cr also includes one resonance, 
etc. Although in principle one should consider all the terms compatible with symmetry, most 
of the large-iV^ phenomenological calculations consider operators with the minimal number of 
derivatives [3S]. This is usually justified through the argument that higher derivative operators 
tend to violate the asymptotic high energy QCD behaviour [HI [39]. Likewise, its has been proven 
in several cases that higher derivative resonance operators can be removed from the hadronic 
action through meson field redefinitions in the generating functional [301 ESI ES EH US HE] • In 
the present article, the leading lagrangian will only contain operators at most 0(p 2 ), with the 
external sources counted as t> M , a M ~ 0{p) and x ~ @(p 2 ) [IS1 HZ]- 

The Lagrangian with only Goldstones has the same form as in xPT but the coupling constants 
are different. In %PT we have the leading order Lagrangian 



C 



(7) 



In R^T beyond leading order the constants standing in front of the operators ( u^u^ ) and ( x+ ) 
may not be the same as in xPT. Therefore, generally we can write 

£ G = ^(n%) + ^( X+ ) (8) 



where we explicitly distinguish between F and F. These can be split in the way, 

F = F + 5F, F = F + 5F 



(9) 



where at large Nc one has the matching condition F = F = F and, hence, 8F and 8F are NLO 
in 1/Nc- On the contrary to what happens in xPT, where the parameters (F and -Bo) which 
characterize the terms ( u^u^ ) and (x+) do not become renormalized, in R/\T the couplings of 
these two operators are needed to make the physical amplitude finite. For simplicity, we choose 
to keep the definitions of the chiral tensors unchanged and to renormalize instead F and F, as it 
was done in Refs. [331 US] with the notation a.\ = F 2 /4 and = F 2 /4. 

The Goldstone bosons couple to massive U(3) multiplets of the type V(l ), A(l ++ ), S'(0 ++ ) 
and P(0 _+ ). The vector multiplet, for instance, is given by 
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(10) 



where we use the antisymmetric tensor formalism for spin-1 fields to describe the vector and 
axial- vector resonances [SI El [32] ■ 
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The resonance fields R are chosen to transform covariantly under the chiral group as in 
Eq. ([5]) [8]. The free- field kinetic term is given by the operators 

£rr = -\{^R»uVaR m ) + \M 2 R (R,uRn + l(V a RV a R') - l -Ml t {R'R') . (11) 

where R = V, A are vector and axial vector resonances and R' = S,P are scalar and pseudoscalar 
resonances. 

The interaction terms which are linear in the resonance fields can be obtained from the seminal 
work [8]: 

C R = c d (Su» Ull )+c m (Sx + )+id m (Px~) + ^( V ^f+) + |^/-k>1> + ^|0W^>- 

(12) 

For our analysis of the SS — PP correlator, the relevant bilinear terms will be [151 US HE] 

£rr- = zAH^P^K) + Af A ({V^,A^}^) + Af p ({V^,P}^>. (13) 

Only single flavor-trace operators are considered for the construction of the large-iVc la- 
grangian. At tree-level, the octet SS — PP correlator only gets contributions from this kind of 
terms, even at subleading orders in 1/Nc- Operators with two or more traces might appear in 
the vertices of one loop diagrams but, since these multi-trace terms are 1/iVc-suppressed, these 
contributions would go to next-to-next-to-leading order and they will be neglected in the present 
work. 

The previous operators provide an appropriate description of the form factors with two Gold- 
stones or one resonance and one Goldstone in the final state. We will perform our most elaborate 
analysis with the lagrangian Lq + £r + Crr>, with at most two resonance fields. As we will see in 
next sections, the R%T description will progressively approach the actual QCD amplitude as more 
and more complicated operators are added. However, although we expect the contributions from 
the operators with three resonance fields to the LECs to be negligible at our level of accuracy, a 
further refinement is eventually possible by considering these operators Crrir» . 

2.2 Subleading Lagrangian 

At the loop level, one needs to introduce new subleading operators in order to cancel the ultra- 
violet divergences, to renormalize R%T and to make the amplitudes finite. As the leading order 
lagrangian operators are 0(p 2 ), the naive dimensional analysis tells us that at one loop one ex- 
pects to find 0(p A ) ultraviolet divergences, requiring the introduction of NLO counter-terms with 
a higher number of derivatives. 

The new operators with just Goldstone bosons required at NLO are, for the SS — PP correlator 
under consideration, 

^ l b° = Y^ + x^) + ^n(x-(vX-^-))-^ 2 ((vX-^-) 2 ) + f (xl-x 2 -)- 

(14) 
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Though we use the same structure of terms as in xPT, the R^T couplings Li are not the same 
as the chiral LECs Lj. The Lj will contribute at low energies to 0(p 4 ) chiral couplings Lj. The 
latter are dominantly saturated by resonances exchanges, so L, are considered to be suppressed 
and sub leading in the 1/Nq expansion. 

In order to make the resonance propagator finite, one needs to renormalize the mass and wave 

functions (M^ } 2 = M r R 2 + SM R , R {B) = Z\R T ) and to introduce at NLO in 1/N C the kinetic 
operator 

C%? = f^V 4 P>, (15) 

with R = S,P. No terms with vector or axial- vectors are needed for the present NLO analysis of 
the SS — PP correlator. 

Likewise, the renormalization of the vertex functions s(x) — ► S and p(x) — > P at NLO in 1/Nq 
will require of the linear terms, 

£™ = Af 8 (5V 2 X+ ) + ^f 3 (PV 2 X -). (16) 

At NLO in 1/Nq, all these subleading counter-terms can only contribute through tree-level 
diagrams. 

2.3 Equations of motion and redundant operators 

The equations of motion (EOM) of the leading lagrangian are given by [46, 33J, 

V% = + + ... (17) 

V 2 S = -M 2 S S + c d u^ + c mX+ + - (18) 
V 2 P = -M 2 P P + td mX ~ + - (19) 

where the dots stand for terms with vector or axial- vector resonances or sources, two-meson fields 
or with one scalar-pseudoscalar external source and one meson field. 

Since most of the subleading resonance operators are proportional to the EOM, it is possible 
to simplify our new NLO resonance operators by means of appropriate meson field redefinitions,: 



C NLO^ C N L 0,e ff = ^ { SS) + C ^( X l)-C m X S M 2 s (S X+ ) 

+^^(PP)-^(x 2 -)-id m X P MUPX-) + - (20) 

where the dots stand for operators that do not contribute to the SS — PP correlator at NLO. 
After the field redefinition the resonance operators and £r LO disappear and the surviving 

terms in the R/\T lagrangian carry in front the effective combinations, 



^8 + -jC 2 m X s — -d 2 m X P + c m Af 8 — rf m Af 3 , 



H e 2 ff = H 2 + c 2 m X s + d 2 m Xp + 2c m X s 18 + 2d m X^ 
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(M 2 s yff = Ml-XsMl 

(M 2 P ) eff = M 2 P -X P M P , 

c e d f = c m -c m X s M 2 -M 2 X s l8 , 

d e £ = d m -d m X P M 2 -M 2 P X( 3 . (21) 

The Lu and L 12 operators do not contribute to terms which can be relevant to our amplitude up 
to NLO and we will see that they are not present in the final result. 



3 Chiral octet SS — PP correlator 

In the case of S"?7(3)-octet quark bilinears, the two-point Green function SS — PP is defined as 



n 



ab 

s-py 



(p) = i d 4 xe ip - x (0\T[S a {x)S\0) - P a (x)P b {0)}\0) = 5 ab U{p 2 ) 



(22) 



1,...; 



with S a = q^=q and P a = iq^'J^q, being A a the Gellmann matrices (< 

In the chiral limit, assumed all along the article, the low-energy expansion of the octet corre- 
lators is determined by %PT in the form [5], 



n( P 2 ; 



2F 2 
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(23) 
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F 2 



32^ 8 ( Mx ) - 



1 38 
7T 2 



1 - In 
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/4 



+ O (N° c ) 



+ 0{p 4 



where in T 8 = 5/48 [3/16] and Tf 8 = -5L 5 /6 [-3L 5 /2] in S77(3)-xPT [E7(3)-xPT]. Notice that 
in xPT the correlator is exactly independent of the renormalization scale /i x , being its choice 
completely arbitrary. 

In the resonance region, one obtains at leading order in 1/Nc, 



9 R 2 p2 



,1 
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m,i 



ML-P 2 



ML-p 2 



(24) 



where one sums over the different resonance multiplets. The subscript $ in M^, c m ^ and d m ^ 
refers to the coupling of the z-th resonance multiplet of the corresponding kind. The requirement 

of the high energy OPE behaviour II(p 2 ) P ~°° 1/p 6 produces the short- distance conditions [39] 
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m.i 



d 2 ■) = — 



J2 c2 m,i M h- d2 m,i M h = - 



(25) 



In the single resonance approximation (SRA), it is then possible to express c m and d m in terms of 
F and resonance masses, 



M 2 P 



8 M 2 -M 2 



d 2 = — 



2 The tiny dimension four condensate -i? (O 



SS-PP 

(4) ) 



8 M 2 -M 2 ' 
-Y2-kcxsF a will be neglected in this work [39l |49] , 



(26) 
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At low energies, we can match the large- No expression (|24|) with the expression ( 123]) . ob- 
taining the LO prediction for the low energy coupling constants L 8 and C 38 , 



c 2 d 2 F 2 



+ 772 , (27) 



2M| 2M 2 16 \M 2 M 2 
r <4^ 2 F± ( Mf, M%\ 

38 2M| 2MJ, 16M 2 M| V Af| Mf, ) K ] 

For the inputs M s = M P /^2 ~ 1 GeV, one obtains L 8 w 0.7 ■ 10~ 3 , C 38 « 7 ■ 10~ 6 for M s = 
1 GeV. However, one does not know to what renormalization scale /i x these numerical predictions 
correspond. In order to pin down this /i-dependence, one must carry the calculation up to the 
loop level. 



4 One-loop computation in resonance chiral theory 

We follow the renormalization procedure presented in [30]. In general, we will use dimensional 
regularization and the MS — 1 subtraction scheme, usually employed in x?T calculations [21 E]. 
This means we will absorb in the coupling counter-terms the ultraviolet divergent piece from the 
loops, counter-terms, 

2 

+ 1e — ln47r — 1 



Aoo(aO = ^ 



d-4 



d-A 



(29) 



Still, the Goldstone propagator and the Goldstone decay amplitudes will be renormalized in 
the on-shell scheme, as it is done in xPT, in order to ease the low-energy matching of R%T and 
%PT at 0{p 2 ). Everything else will be renormalized in this section in MS — 1. For simplicity, 
we will denote this set of schemes as MS from now on. Afterwards, we will^ study alternative 
renormalization schemes for the R%T couplings and their relation with the MS parameters. 

In this section, together with the general structure of the amplitudes, we will provide in this 
Section just the explicit results for the case when the lagrangian contains the operators Cq + Cr 
with at most one resonance field, derived by Ecker et al. [8|. The contributions from operators 
£rr' with two resonance fields are provided separately later in Sec. [71 For clarity, we provide the 
individual contributions from each absorptive cut (e.g. tttt, Vtt...). The precise definitions for the 
corresponding Feynman integrals are given in Appendix O 

4.1 Goldstone boson renormalizations 
4.1.1 Goldstone self-energy 

The general form of the renormalized Goldstone propagator is given by 

l \ = —j^P p2 ^(P ) > (30) 



with Z<j, the wave function renormalization of the bare Goldstone field, <jr- ' = Z£cjf. In order to 
make the propagator finite, one needs to perform the shifts 

Z^=1 + 5Z^ F = F + 5F, L 12 = L\ 2 + 5L l2 , (31) 
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A 



o- 



Figure 1: Contributions to the Goldstone boson self-energy. The single line represents the Goldstone 
boson while the double line represents the resonance. The type of resonance is written above it. 



where SZ^ and SF are NLO in 1/Nc- The NLO coupling L 12 is split into a finite renormalized 
part L\ 2 and an infinite counter-term 5L 12 . 

Considering the on-shell renormalization scheme for the Goldstone propagator, i.e. such that 



p- 



0(p A ), leads to the renormalization condition 



25F 



+ SZ^ - E» 



(32) 



with E>) - ^ 



dp 2 



. The 0(p 4 ) ultraviolet divergence in E^ is absorbed into <5L 12 in the MS 



scheme. The renormalized Goldstone propagator is then provided by 



p 



W) 



with its perturbative expansion, 

a: = 



P 



+ 



P 



4L[ 2 p 4 



+ K(p 2 



+ 



(33) 



(34) 



where the dots stand for the next-to-next-to-leading order corrections (NNLO) and 
E£(p 2 ) = ^(p 2 ) -p 2 E;(0) - S^ 2 )| Aoc0(p4) behaving like 0(p 4 ) when p 2 - 0. 

If one considers just the contributions Cr from interactions linear in the resonance fields [8], 
the one loop Goldstone self-energy E^ is given by the diagrams shown in Fig. HJ A priori, tadpole 
diagrams might appear, either with a Goldstone or a resonance running within the loop. However, 
they happen to be zero in the chiral limit. All this yields the renormalizations and the renormalized 
self-energy, 



25F 
~F~ 

5L 12 



+ SZ 6 + 



8F 4 7T 2 



9GyMy 



Aoo + In ■ 



V 2 



+ 6 



„ ., / , Ml 1 

3r-, .V s - ( Aoo+ln-f -- 



3(2c 2 + G 2 



64vr 2 F 2 



0. 



(35) 



with 



<j>(x) = 



2 m i 



\S<f> 



v<t> 



3cjP 
8vr 2 F 4 

ZG 2 v p A 

16tt 2 F 4 



, M s 
ln^ + 



In ■ 



M, 



v 



P z 



+ 



P 

Ml 
P 2 



My 



1 - 



x 



ln(l — x) — 



(x - 2) (x - |) 



x- 



(36) 



(37) 
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- - A 

L 1V L 12 , 8 F, 8 Z A 



go- 



Figure 2: Contribution to the vertex p<j). The crossed circle stands for a pseudo-scalar density insertion. 



4.1.2 Vertex p<$ 

The vertex function has the form 



<f> p<t> (p 2 ) = V2 



ZJF 2 B 4V2B< 



oP 



(L n + L 12 ) + <& p(j> (p 



2\1£ 



(38) 



where Q p( p(p 2 ) u represents the one-particle-irreducible (1PI) contribution from meson loops. 

Notice that it is convenient to choose the renormalization scheme for SF such that the on-shell 
decay amplitude coincides with the pion decay constant, which by construction we denote as F. 
Thus, for the renormalizations 



one has 



F = F + 8F, 



— + 2^ 



Ln = L u (/i) + <5Ln(/i) 



V2B F 



0. 



(39) 



(40) 



and the counter-term 5Lu(/j,) is chosen to cancel the 0(p 2 ) divergent terms in & pc f>{p 2 ) u in the 
MS-scheme. The renormalized vertex function is then equal to 



V) = ^ b of\ i - ^ _ + ^%,(pY 



F 2 V2B Q F 



(41) 



with $^(p 2 ) 1£ being 0(p 2 ) when p 2 — > 0. 

In the case with only Cr interactions, linear in the resonance fields [5], one has the diagrams 
shown in Fig. [2l These lead to the renormalizations and renormalized one-loop contributions, 



25F 1 



SZ, 



F ' 2 
SLu(fi) + 6L 12 (n) 



0. 



16tt 2 F 2 



-^00 > 



(42) 
(43) 



with 



- $ r Jp 2 ) u \ 

V2B F pAP> 



ip(x) 



Stb 



3c d c m p 2 
4tt 2 F 4 



1 , M s , f 



1 

X 



ln(l 



x) 



(44) 



(45) 
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Figure 3: Contribution to the vertex pcj>. The crossed circle stands for a axial-vector current insertion. 



4.1.3 Vertex a</> 

Although it is not required for the correlator calculation in this article, we will compute the a<p 
vertex function for sake of completeness. From previous calculations we obtained two equations 
for three unknown objects SF, 5F and SZ^. The third equation can be found by analyzing the 
a M — > (j) vertex, which, abusing of the notation, has the form 

<Mpf = 3W • V* (46) 

where 

*^V2M-i^ + ,Mpy. (47) 

As it happened before with SF, it is convenient to choose for 5F (as we did here) the scheme 
that recovers the pion decay constant F when the decay amplitude is set on-shell {p 2 — > 0): 

+ \m* + = °- ( 48 ) 

The coupling £Li 2 (/i) is chosen to cancel the 0(p 2 ) UV divergent term in <& a <t>{p 2 ) 11 i n the MS 
scheme. 

When only Cr interactions are taken into account [8], the diagrams shown in Fig. [3] yield the 
renormalizations 



2SF \ grw 1 
+ -SZ 6 + 



F 2 8vr 2 F 4 



VGlMlf, , M 2 1\ o2 „, 2 A , M 2 1 
' 1 ' Aoo + In + - - 3c 2 ,M 2 + In 



(49) 



SLM = -?M±S) Aoo . (50) 

In this case, it is possible to see explicitly that the renormalization for Lyi is in an agreement with 
its former result from the Goldstone propagator. 

4.1.4 Renormalization of F, F and SZ^ 

Comparing the three equations for SF, SF and SZ^, one is finally able to extract each of them 
separately: 

6Z+ = 2E;(0)" + ^<M0)" , 
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5 M s , X s , Z s 



Figure 4: Contributions to the scalar resonance self-energy 



if 



V2F 



<Mo) 



if 



(51) 



Thus, in the case when only interactions Cr, linear in the resonance fields, are considered [5], one 
gets 5Z& = 0, 5F = and 



5F 



16tt 2 F 4 



•«W Aoo+ln 4 i 



3c 2 M| ( Aoo + In ^# - ~ 



Ml I 

2 2 



(52) 



This confirms the results from Ref. [33J, where F was renormalized but F was not. On the other 
hand, the renormalizations of F and F were not considered in Ref. [30] and, consequently, a 
nonzero 5Z,h was found. 



4.2 Scalar resonance renormalization 
4.2.1 Scalar resonance self-energy 

The renormalized propagator has the form 

iA~ s l = Z s (p 2 - M 2 S ) + X s p* - S 5 (p 2 ) 



(53) 



where we have performed the scalar resonance wave-function renormalization = ZgS r . In 
order to cancel the divergent terms of the one-loop self-energy £s(p 2 ), we make the shifts 



M 2 s = M r s 2 + 5M 2 S , Z s = 1 
The renormalized propagator is then given by, 



SZ, 



X s = X s (fi) + SX s {^. 



^s 1 

with its perturbative expansion, 

i 



P 



Ms + X's(fi)p - TT s {p 2 



A, 



+ 



p 2 -M r s 2 (p 2 - M r s 2 ) 2 



- X s {^ + E£(p 2 ) + 



(54) 
(55) 

(56) 



In the case where only the Lr interactions are considered, one obtains 
5M S = 0, 5Z S = 0, 6X s (ji) 



3ci 



16vr 2 F 4 



2^,4 



3cjp 
16tt 2 F 4 



l - In 



~P 

V 2 



(57) 
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Figure 5: Contributions to the vertex sS. The crossed circle stands for a scalar density insertion. 
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Figure 6: Contribution to the pseudoscalar resonance self-energy 
4.2.2 Vertex sS 

The vertex function s(x) — ► 5 has the form 

$ s5 (p 2 ) = -4P | z\c m - Af 8 p 2 - ^^(p?)" J . (58) 

The renormalizations of the scalar wave-function Z$ = 1 + SZg, the LO constant c m = c r m {n) + 
<^c m (//) and the NLO coupling Af 8 = Af 8 (/i) + <5Af 8 (/i) make the amplitude finite: 

$ sS (p 2 ) = -AB | C - Af 8 (/V - ^^(p 2 )" } • (59) 

In the case with only Cr interactions [Hj, we had SZs = 0. The cancelation of divergences in 
the MS scheme leads to the shift and the renormalized one-loop contributions, 

*c m = 0, 5 x s M = -^ ¥ - 2 X O0 (60) 

l <s(P 2 ) U U = ^fl-ln^V (61) 



AB sb ^ ' lw 64tt 2 F 2 V /u 2 

4.3 Pseudo-scalar resonance renormalization 
4.3.1 Pseudoscalar resonance self-energy 

The renormalized pseudoscalar propagator has the form 

zAp 1 = Z P (p 2 - Mp) + X P p 4 - S P (p 2 ), (62) 

i 

with p( B ) = ZpP r . The cancelation of the Aoo UV divergent terms in the one-loop self-energy 
Ep(p 2 ) needs the shifts 

Mp = M r P 2 + 5M P} Zp = l + 5Z P , Xp = X P (n) + 5Xp(Li), (63) 
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Figure 7: Contribution to the renormalization of vertex pP 

leading to the renormalized propagator, 

zAp 1 = p 2 - M P 2 + Xp(^p 4 - S P (p 2 ) , (64) 
and its perturbative expansion, 

A ' = w^w + w^w? { - XUp)pl + E ^ 2 > } + • • • < 65 > 

In the case where only the Cp interactions are considered [8], there is no one- loop diagrams 
contributing and, therefore, SZp = 5Mp = 5Xp = 0. 

4.3.2 Vertex pP 

The vertex function p(x) — > P has the form 

% P (p 2 ) = -4B Q | 4 dm - Af 3 p 2 - ^^(P 2 ) 1 ' } • (66) 

The renormalizations of the scalar wave-function Zp = 1 + 5Zp, the LO constant d m = d r m {n) + 
Sd m (fj,) and the NLO coupling Xf 3 = Af 3 (/i) + 5Xf 3 (fi) make the amplitude finite: 

% P (p 2 ) = -AB Q | C - Af 3 (//)p 2 - -^%p(p 2 ) U } • (67) 

In the case with only C R interactions [SJ, £Zp = and there is no loop diagram contributing 
to this vertex, so we have 5d m = 5Af 3 = and the renormalized vertex function results 

%p(p 2 ) = -4£o R - Af 3 ( /U )p 2 } . (68) 

4.4 1PI contributions 
4.4.1 1PI diagram ss 

Now, we analyze 1PI diagrams that appear in the ss-correlator: 

IOT) = 16S Q % + 8B 2 H 2 + ICV) 1 *- (69) 

The shifts L 8 = L%(fj,)+8L$(n) and if 2 = H^ifj) + 5H 2 (fi) render the amplitude finite by canceling 
the UV divergences in the MS'-scheme, which becomes 

nL p V) = i6£ 2 W + 8B 2 H r M + n^V)* (70) 
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Figure 8: Contribution to the 1PI vertex ss 
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Figure 9: Contribution to the 1PI vertex pp 



In the case with only interactions Cr linear in the resonance fields [8], the 1PI diagrams 
contributing to the S'S'-correlator are shown in Fig. [HJ Thus, one gets for the shifts and the 
renormalized amplitude the expressions 



25LM + 5HM 



3(F 2 + 16dl 
128vr 2 F 2 



-Xr 



(71) 



n lPI,r/ 2\U\ 
^ss VP ) \H> 

VV ss VP ) \P<t> 



B, 



167T 2 

2 oa m 
° 7T 2 F 2 



In 



In 



-P 

Ml 



Mj 

~7 



In 1 



P 

Ml 



(72) 



4.4.2 1PI diagram pp 

Similarly, for pp-amplitude one has the structure, 



CV) 



16£? 2 L 8 - 16B 2 L U 



8B 2 L l2 + 8B 2 H 2 + TlHW 



(73) 



The UV divergences are absorbed through the renormalization of L 8 , L 11; L 12 and H± 2 , rendering 
the amplitude finite: 

n^V) = -165 2 Z^(/i) - lQB 2 U n ^) - 8B*L[M + 8B 2 H r M + U^(p 2 ) u . (74) 

In the case where only the contributions from Lr operators are considered [8], the divergences 
are absorbed by the shift 



26 LM + 26L u (n) + 6L u (n) - SH 2 



V 2 

87T 2 F 2 °°' 



leaving the finite one-loop contribution, 



lPLr/ 2\U\ 



B 



2 OL m 

^2 T72 
7T r 



In— £ 

/i 2 



P 



In 1 



P 

M 2 



(75) 



(76) 
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4.5 Correlator at NLO 

At NLO we can write the general 1PI decomposition of the SS — PP correlator in terms of 
renormalized correlators and vertex functions, 

U ss _ pp (p 2 ) = tA s (p 2 ) (<MP 2 )} 2 " iA P {p 2 ) {%p{p 2 )} 2 - i^{p 2 ) {%^P 2 )} 2 

+ - I4 P V), (77) 

where we made use of the relation between the vertex functions for incoming and outgoing mesons, 
$ s5 = $ 5s , <S> pP = $ Pp , <$> p4> = $ # . 

If one now uses the previous perturbative calculation, the SS — PP octet correlator takes up 
to NLO in 1/N C the form, 

2 v „4' 



-Lu(p 2 ) = 2 1 J 16c4 - 32c m Af 8 p 2 + 16c ™ XsP ) 
#o M 2 s -p 2 \ m 1SF Ml-p 2 J 



(M 2 -p 2 ) 2 J M 2 -p 2 B 
1 (l6d 2 m -32d m Xr 3 p 2 + 1 -^4 



Mp - p z V Mp — p 



r 2 „2n2 2j pIP ) + s>f 2 2 d ^pP\P ) 



(Mp-p 2 ) 2 ^ ' Mp-p 2 B 

+32L 8 + 16L n + 8L 12 + U r ss _ pp (p 2 ) u . (78) 

The couplings shown here (and from now on) are the renormalized ones even if the superscript 
"r" is not explicitly present. The first two lines are the contribution from the scalar exchanges. 
The third and fourth ones come from the pseudoscalar resonance exchanges, whereas the fifth 
one is produced by the Goldstone exchanges. The last line is given by the 1PI diagrams in the 
SS — PP correlator. 

Notice that the correlator results independent of L\\ and L 12 due to the cancelation between 
the Goldstone exchanges and the 1PI terms in f!78|) . Likewise, it is possible to check that the 
correlator only depends on the effective combinations c^, d^, M| ff , MJ? L| ff from Eq. (I2TT) : 

1 1 fir eff 2 1 fir 2 Rr 1 

B 2 yF! Mf 2 -p 2 {M 2 s -p 2 ) 2 SKF J Mg — p Bq sSKy j 

16^m . 16(i m l2\\l , &d m 1 /2\U 



I _ m yr ( 2\W , — '» - for / 2 

reff 2 2 + / ji,f2 J2\2^P\P ) + A# 2 ^2 R ^pPlP 



M P tt2 -p 2 (Mp-p 2 ) 2 ™' M 2 -p 2 B 
2F 2 2F 2 2 ,» 2F v2_ r , 2 \i/ 

+32Lf 2 + nL_> 2 )^. (79) 
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The couplings Xs, Xp, Af 8 and Af 3 disappear from our NLO calculation and c m , d m , Mr and 
L 8 are replaced everywhere by c^f, d^, M R S and Lf 1 . The replacement in the subleading terms 
leaves the expression unaltered up to the order in 1/Nq considered in our computation. 

This elimination of the renormalized couplings X s , X P , Af 8 and Xf 3 can be understood in 
an equivalent way by means of the EOM of the theory and the meson field redefinitions. The 
effective couplings that are left in front of the operators after the meson field transformations 
coincide exactly with the combinations that determine the correlator up to NLO. 

In the subleading terms in Eq. (1791) . a priori one can use indistinct the original couplings, 
e.g. c m , or the effective ones, this is, c^f, as the difference goes to NNLO. However, for sake 
of consistence, one should always consider the same renormalized coupling everywhere in the 
amplitude. Hence, after performing the field redefinition that removes Xs, Xp, Af 8 and Af 3 , all 
the remaining couplings appearing in H(p 2 ) are the effective ones. From now on, we will consider 
that the R/\T action has been simplified through meson field redefinitions in the previous way and 
the superscript "eff" will be implicitly assumed in the couplings in order to make the notation 
simpler. 



5 High energy constraints 

The NLO expression for the correlator contains plenty of resonance parameters that are not fully 
well known. A typical procedure to improve the determination of these couplings is the use of the 
short-distance conditions j9]. 

The operator product expansion tells us that the SS — PP correlator vanishes like 1/p 4 for 
the large Euclidean momentum. Indeed, due to the smallness of its dimension-four condensate 
(^{Of~ pp ) ~ UnasF 4 ~ 3 • 1(T 4 GeV 4 [39]), it is a good approximation to consider that it 

vanishes like 1/p 6 when p 2 — > — oo [39, 49J. 

The R%T correlator does not follow this short-distance behaviour for arbitrary values of its 
couplings. This imposes severe constraints on the coefficients of the high-energy expansion of our 
NLO correlator, 

4n(p 2 ) = E w + ^^f) ■ ( 8 °) 

^0 n=0,l,2„. ) V H- J 

The proper OPE short- distance behaviour is therefore recovered by demanding [17] 

af = = , for k = 0, 2, 4 . (81) 

At large Nc, there are no logarithmic terms (cxf* = 0) and for the remaining coefficients 
one has aty = (no L 8 or higher local couplings at large iV^) and the two Weinberg sum-rules 
(WSR) [3S], 

c4 p) = 2F 2 + lQd 2 m - 16(4 = 0, 

= \U 2 m M 2 P - lQc 2 m M 2 s = 0. (82) 

At NLO, in the case when the interactions only contain operators Lr with at most one reso- 
nance field [8], the high-energy expansion log-term coefficients result 

= 8c 2 m -8d 2 m -4c d c m + 2c 2 d + G 2 v - 8 -^-^, 
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-Olr 



8d 2 m M 2 - 8c 2 m M 2 s - lbM f 2 d m + 20c d c m Ml - 6c 2 d M 2 s - 3G 2 V M 2 , 



F 2 

^ ' 2 W - 24C ^ C - M| - 4c d c m Ml + 6c 2 Ml + 3G V M V , (83) 



-a 



i 



F 2 

and the high-energy coefficients 4 are given by 

= -a { n ] + 32 Z 8 



a 2 W = 2^ + 16(4-16(4 + 

af = lQd 2 m M 2 - 16e 2 m M 2 + B(n), (84) 



with the NLO corrections 



tt 2 F 2 V /" 2 7 ^ 2 V J k 2 F 4 

6c d c m Ml f^Ml 1\ 9c 2 M 2 /M 2 A 9G\M V / M 2 , 1 



7T Z 



F 2 V 4y 4vr 2 F 2 V H 2 2 J 8vr 2 F 2 V A* 2 2 



3(e.Mj> , Qc^Mj , 3ciMj 6c d c m M| Q^M* / M 2 1 



W - 2 F% 2 F% 2 2F% 2 tt 2 F 2 4vr 2 F 2 \ /x 2 2 
, 3Q Cm M| M 2 9G V M V ( M 2 1 



tt 2 F 2 /i 2 8vr 2 F 2 \ /j 2 2 



^-T--^f -- . (85) 



The large-iVvj WSR (|25l) gain the subleading contributions in 1/Nc, yielding for af = oi 
the solution! jS], 



4 



_ F^ M| / A(/x) _ gfr) \ 2 _ F^ M| / Afc) 



8 Mp — M| V 2F 2 2F 2 M 2 J' m 8 M 2 -M 2 \ 2F 2 2F 2 Mj 



(86) 



The couplings Ms, Mp, c m and d m may also depend on /i. Nonetheless, unless necessary, this 
dependence will not be explicitly shown. Also, as stated at the end of the previous section, 
one must keep in mind that these are the results after the meson field redefinition that removes 
the redundant couplings Xs, Xp, Af 8 , Af 3 , so the surviving couplings carry the superscript "eff" 
implicit. 

The «o = constraint implies that L 8 = also at NLO in 1/Nq (for any renormalization 
scale n). We will see that for all the possible interactions considered in this paper, now here and 
later on, there is the same constraint = — aig +32L 8 and, therefore, in general we find L 8 = 0. 
The constants af\ A(fi) and F(/i) only arise at NLO or higher. Hence, when they are used for 
the computation of the correlator up to NLO, one can indistinctly use for their calculation either 
renormalized couplings or their large- Nc values, as the difference goes to NNLO. 

Although these expressions will be used later in other renormalization schemes, A(fi) and 
B(n) will always refer to their former definitions in the MS scheme, like, for instance, the results 
provided in Eq. fl85l) . 



3 The notation A(jjl) = 2F 2 5 { ^ LCn 2F 2 M^ ] LO = B(p) was used in Ref. [31] 
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5.1 Alternative renormalization schemes 

During the renormalization procedure we considered the MS'-subtraction-scheme for all the res- 
onance couplings. However, in some situations one may get large contributions from A(fi) and 
B(n). The NLO prediction for c m and d m derived from Eq. (IBS"]) may then become very different 

from the large-Afc WSR determinations c 2 m = ^ m^-m* > ^li = x m^-m^ ■ 

A way out to minimize possible large radiative corrections to the WSR is the choice of con- 
venient renormalization schemes for couplings (c m and d m ) and masses (Ms and Mp). In the 
renormalization procedure we originally chose to cancel the Aqo from the one-loop diagrams, but 
we could have chosen to cancel the Aoo term plus an arbitrary sub leading constant. This change 
makes that instead of having in the amplitudes the renormalized coupling AIL^ in the first scheme, 
one now has the renormalized coupling in the second scheme plus a constant, A^ 2 + C #1_>#2 . Thus, 

effectively one can account for a change from the MS'-subtraction-scheme (with renormalized 
couplings k = c m , d m , Mf , Mp) to another (with parameters k = c m , d m , Mj, Mp) through the 
shifts, 

k = k + Ak. (87) 

The difference Ak will be, of course, subleading in the 1/Nc counting with respect to k and k. 
This will affect the parts of the calculation where these couplings contribute at LO in 1/N C - In 
the contributions that start at NLO (e.g. A(n) and B(fi)), the variations due to Ak go to 
NNLO and they are therefore neglected. If one applies this change of scheme to Eq. (186]) . one gets 
for the NLO extension of the WSR, 

o> } = 2F 2 + 16d 2 m -16c 2 m + (32d m Ad m -32c m Ac m + A(^ , 

a ( l ] = 16 d\Ml- 16 c 2 m M 2 s (88) 
+ (32 M 2 P d m Ad m + 16 d 2 m AM 2 P - 32 M 2 s c m Ac m - 16 c 2 m AM% + £(//)) . 

The terms within the brackets, (•••), would be the finite contributions from the one-loop diagrams 
in the new scheme. 



5.1.1 Pole mass scheme for M s and M P 

In addition to the MS-scheme for the scalar and pseudo-scalar masses (AM^ = 0), we will also 
study the pole-mass scheme. The problem with the MS mass is the difficulty to give a direct 
physical meaning to the /^-dependent mass M^(/i), specially when more and more operators are 
added to the R%T action. On the other hand, the resonance pole mass is a universal property 
which does not rely on any particular lagrangian realization. Thus, instead of considering the \x- 
dependent renormalized masses Mp(fi), we will switch to the renormalization scale independent 
pole masses M R = Mp lc , defined by the pole positions (M^° lc — ir^° lc /2) 2 of the renormalized 
propagators. Up to NLO in 1/Nc, one has 

MT 2 = M 2 R + Re^ R (M 2 R ), M^rT = -l^ R (M 2 R ) , (89) 

and therefore, 

AM 2 R = M 2 R -M 2 R = -ReZ R (M 2 R ), (90) 
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Since AM^ is NLO in 1/Nc, the difference between using the M R (MS'-subtraction-scheme) 
within E(M|) or its value Mr in another scheme goes to NNLO. Therefore, it is negligible at the 
perturbative order we are working at. 

If only interactions Cr given by operators linear in the resonance fields are taken into ac- 
count [8], one has for the pole scheme 



16vr 2 F 4 



1-ln— £ 



AMp 



0. 



where only the two-Goldstone loop £s(p 2 )|< 
interactions are taken into account [81. 



contributes to AMj and Sp(p 2 



(91) 

if only the Cr 



5.1.2 WSR-scheme for c m and d m 

Since the value of the spin-0 parameters is very poorly known at the experimental level, one 
finds important uncertainties and variations in the determination of c m and d m through the NLO 
sum-rules (1861) . The choice of a shift that minimizes the finite part of the loop contributions is 
not straight-forward. For instance, within the MS-subtraction-scheme itself, it is not easy to 
find a value of \x that minimizes both A(fi) and B(fi) at once unless the resonance couplings are 
appropriately fine-tuned. This makes the short-distance matching rather cumbersome and the 
extraction of the necessary resonance parameters problematic. 

Alternatively, the selection of a shift An that exactly cancels the one-loop contributions to 
Eq. flHBT) (provided in the MS'-scheme by the constants A(/a) and B(fi)) seems to be a better 
option. This converts Eq. (|86|) into 

a (p) = 2F 2 + lg ^ = o , a ( t = 16 d 2 m M 2 P - 16 c 2 m M 2 s = , (92) 

with the solutions 

£2 _ F 2 M 2 g _ F 2 Ml 

° m ~ T M 2 -Mf m ~ ~8 ' M 2 -M 2 (93) 



Though this has the same structure as the LO prediction (fStjj) of the large-iV^ WSR in Eq. 
the couplings appearing here are the renormalized ones. Nonetheless, this result ensures that the 
difference between c m /F and d m /F at NLO and their large- Nc limits remains small provided 
Mr w Mp c ^°°. In order to achieve this minimization, the shifts An must be tuned in such a way 
that they obey 

32d m Ad m — 32c m Ac m + A(/i) = , 
32Mpd m Ad m + lGd^AMp - 32M|c m Ac OT - 16c^AM| + B{ji) = 0. (94) 

If one fixes AMr (for instance, through the pole scheme) the solutions for Ac m and Ad m are then 
given by 

MpAjp) - B(n) + 16 c 2 m AM 2 - 16 d 2 m AM 2 
m m M 2 - Ml 

32 J Ad MlAjp) - Bjfj) + 16 c 2 m AM 2 - 16 djAM 2 

m m M 2 -M 2 ' { ' 



23 



In the change of scheme we will make the replacement k = k + An in the tree-level LO 
diagrams, whereas in the subleading contributions we will just consider the difference 

goes to NNLO in 1/Nc- Thus, we will end up with a matrix element expressed in terms of just 
renormalized couplings in the new scheme (k). We will denote the c m and d m renormalization 
scheme prescribed by Eq. fl95l) as WSR-scheme. 

6 Low-energy expansion 
6.1 MS— subtraction scheme 

At low energies, the expansion of our one-loop R%T correlator yields the structure, 
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where in the U(3) case we obtain G 8 = jq 

using the LO matching relation L5 = CdC m /Mg [8]. The logarithm from the tttt loop in R%T 
has been singled out in the ln(— q 2 ) terms. These R%T logarithms exactly reproduce those in 
the low-energy x?T expression (1231) . ensuring the possibility of matching both theories [47J. The 
one-loop contributions from the remaining channels generate only polynomial terms at this chiral 
order and they are provided here by £l 8 and £c 38 , defined within the MS'-renormalization-scheme. 
The predictions for the low-energy constants at NLO in 1/Nc then turn out to be 
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The dependence of the terms on the right-hand side of the equations (r.h.s.) on the R%T renor- 
malization scale /i have been left partially implicit. Only the last term shows \x explicitly. It 
comes from the two-Goldstone loop in R%T (Eq. (1961) ) and matches exactly the log from the 
two-Goldstone loop in %PT (Eq. (I23l . with the chiral renormalization scale /i x ), producing the 
ln(/i 2 //x 2 ) term. This ensures the right low-energy running with \x x for the ;^PT low-energy con- 
stants [17]. On the other hand, the r.h.s. is independent of the R%T scale \x at the given order 
in 1/Nc- There can still be some residual \i dependence at NNLO, which would allow the use 
of renormalization group technics in order to improve the perturbative expansion and to remove 
possible large radiative corrections [38] . Nonetheless, this is beyond the scope of this article, where 
we will take the usual prescription fj, — \i x (30l 061 Wf\ . 

If we use the c m and d m predictions from the high-energy OPE constraints in Eq. ( |86l) . the 
low-energy predictions result [31] 
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In the case, where we only have interactions Cr in the lagrangian, linear in the resonance 
fields [H], the low-energy contributions from the one- loop diagrams are given by 
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The results ( |98|) correspond to the predictions for the U{3) chiral perturbation theory couplings, 
where the T]i is identified as the ninth chiral Goldstone. In order to recover the traditional SU(3) 
couplings one needs to make use of the matching equations [31, 50J, 
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These outcomes will be used later in the alternative renormalization schemes andjihe constants 
CisO")) £c7 38 (/f)) ^4(a*) and will always refer to their former expressions in the MS scheme. 

6.2 Pole masses and WSR— scheme for c m and d m 



In this case, the renormalization scheme of c m and d m is chosen such that the one-loop contributions 



to the NLO relations in Eq. (1561) are exactly canceled, yielding c; 
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The low energy limit of the R%T correlator in the new scheme leads to the LEC determination, 
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AMJ(MJ + 2MJ,) AMJ,(2MJ + Mp) 
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where £l 8 ,c 38 are the same one- loop contributions to the LECs computed before in the MS- 
subtraction scheme. The same applies to A(fi) and -B(/i), which were defined as the one-loop 
contributions to the high-energy expansion coefficients in the MS-scheme. In £ Lg , £c 38 , and 
_B(/i) we will use the couplings and masses in the new scheme (c m , d m , Mr) instead of the original 
ones in the MS-scheme (c m ,d m , Mr), as the difference goes to NNLO in 1/Nc- The constants 
AM| = M R — M\ provide the difference between the mass Mr in the MS'-scheme and its value 
Mr in another scheme. In this paper it will refer in particular to the mass pole, although it accepts 
further generalizations. 

Notice that these expressions are similar to those in the MS'-scheme ( 198|) . up to the AM| 
terms that arise due to the change of mass prescription. The WSR-scheme does not modify the 
low-energy prediction, it just serves to reduce the uncertainties in the NLO Weinberg sum-rules. 

Finally, in order to obtain the traditional SU (3)-%PT LECs, one should use again the matching 
Eq. (PHDD- 



7 Correlator with the extended RxT lagrangian 

Ecker et a/.'s lagrangian [8] has been found to be very successful for the description of amplitudes 
with few-Goldstones (im form-factors, scatterings...). However, it fails to describe processes with 
multi-Goldstones states or with a higher number of resonances. The LO meson lagrangian must 
be then enlarged to improve the description of the new channels. In the case of our observable, 
the relevant operators with two resonance fields are [I5j [32, |33l HE] , 

£rr> = ^\^( [V"P, V^} u v > + Af A ( {WS, A,,} u v ) + Af p ( { WS, P}u,). (102) 

The Af v and Af p terms induce a one-loop mixing between the Goldstone and the pseudoscalar 
resonance. These loops bring ultraviolet divergences which need the presence of the subleading 
counter-terms, 

AC P = d'J PVX )+d'L( (V 2 P)V/ ) , (103) 

to make the amplitude finite. At LO, in the free field case, the meson kinetic terms are assumed 
to be defined in the canonical way, i.e., without mixing between particles. This was indeed the 
case in Ecker et a/.'s lagrangian [8]. In addition, although these P— <f> operators may arise at NLO, 
they happen to be proportional to the EOM. They can be removed from the action through a 
convenient meson field redefinition, leaving for the relevant couplings in our problem the effective 
combinations 

Lg ff — L 8 + -c 2 m X s - -d 2 m X P + c m Af 8 - <i m Af 3 - -d m d'm > 

K\ = H 2 + c 2 m Xs + d 2 m Xp + 2c m Af 8 + 2d m X^ 3 + d m d n m , 

(M 2 s ) eff = M 2 s -X s Ml, 

{M 2 P ) eff = M 2 p-X P Mp, 

c e l f = c m -c m X s M 2 s -M 2 s X s l8 , 



26 



o --o 

PV, SA, SP 



Figure 10: Contribution from Crr> operators to the Goldstone boson self-energy. 



dt S = d m -d m XpM 2 P -M 2 P \^ + ^d' m -^M 2 P dl. 



(104) 



7.1 Meson self-energies 

These operators do not modify the previous loop contributions. However, new channels are now 
open in the different vertex-functions. Thus, the Goldstone self-energy gains the contributions 

(Fig. Eon, 
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in addition to the former V(j) and S<j) cuts from Eq. fl36l) . The functions J(p 2 ,M 2 ,M 2 ) is the 
subtracted two-propagator Feynman integral ( J(0, M 2 , M 2 ) = 0), given in App. O 

The scalar propagators contains now Acf) and P<p cuts (besides the <jxft-one from Eq. flBT|) ) 
(Fig.HU): 
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Figure 11: Contribution from £pp' operators to the scalar and the pseudo-scalar resonance self-energies. 



Ecker et al's lagrangian C R did not modified the pseudoscalar resonance propagator. However, 
the new operators Crw yield (Fig. [Til . 
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The renormalized resonance self-energies provide at this order the pole masses through Eq. (J89 
giving the corresponding shifts AMj, = — Re£p(Mp). 



7.2 P— (f) mixing 

In addition, these operators Af y and Af p also generate a P—<p mixing (Fig. I12p . 
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leading to an extra perturbative contribution to the PP-correlator that has to be added to the 
former ones in (I77p : 
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After a convenient field redefinition d' m and c?^ disappear from Eq. (11091) . being their information 



encoded in df , Lf and iff. 

It is important to remark that at the NLO under consideration, the mixing does not modify the 
pseudoscalar resonance mass renormalization. The Goldstone remains massless -as expected- and 
the resonance pole mass is still provided at this order by AMp = — ReS P (Mp) through Eq. fl89|) . 
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Figure 12: Contribution from Crri operators to the mixing term between the Goldstone and the pseudo- 
scalar resonance. 



Figure 13: One-loop diagrams with Crri operators in the s(x) — ► S and p{x) — ► P vertex functions. 

7.3 New s — > S and p — > P vertex functions 

A new Pip channel is opened in the s — > S vertex function in addition to the 00-cut from Eq. (IBTj) : 



8 Phenomenology 

The R/\T lagrangian developed by Ecker et al. [8], £ = Co + £r, only contained operators 
with at most one resonance field. This approach has been proven to be very successful at the 
phenomenological level for the last two decades |39j . Nevertheless, in the few last years it has 
become clear that the description of more complicated QCD matrix elements (e.g. 3-point Green 
functions [HI [151 EH El EES] ) demands the introduction of operators with more than one resonance 



Since the Mr(//) masses in the MS — scheme are \x dependent, they are difficult to relate with 
the physical masses provided, for instance, by the Particle Data Group (PDG) [51J. This relation is 
even more cumbersome when one adds more general kinds of vertices (e.g. Af p ) within the loops: 
in the MS'-scheme the value of Mr(//) will depend on the content of the theory and its lagrangian. 
Thus, it seems more convenient to use universal properties such as the pole masses, denoted here 
as Mr. The octet of the lightest scalar and the pseudoscalar resonances are then related, to the 
ao(980) and the 7r(1300), and we will consider from now on the inputs Ms = 980 ± 20 MeV and 
M P = 1300 ± 50 MeV [281 EI]. 





not get any contribution from Cr alone: 




field [T5]. 
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The procedure that we will follow in order to extract the LECs with higher and higher accuracy 
is to progressively add more and more physical information to the RxT correlator, starting from 
lower energies. Since the resonance parameters will be used to accommodate the short-distance 
OPE behaviour, in general the two-meson thresholds (Sir, Vtt, Pit...) may not be at the right 
place. Likewise, one may find that individual intermediate two-meson channels have a clearly 
erroneous momentum dependence at high energies (e.g. constant or growing behaviour). 

The introduction of the new operators \\ p , Xf p and Xf A will allow us to improve the momen- 
tum dependence of the Rtt absorptive channels with one resonance and one Goldstone. However, 
since these new couplings will be tuned to implement the short-distance OPE constraints, the 
Rtt channel description may still differ slightly from that provided by the physical values of Xsp, 
Xpv, Xsa, Cd, Cry... Likewise, the two-resonance RR' absorptive cuts will still remain wrongly 
described until operators with three resonance fields are taken into account. Nonetheless, we will 
see that the R%T description progressively approaches the actual QCD amplitude as the hadronic 
action is completed with more and more complicated operators, bringing along a better and better 
description of the lower channels. 



8.1 Phenomenology with Ecker et aUs lagrangian Cq + Cr 

First, we will extract the value of the LECs at large Nq within the single resonance approximation. 
We will use the formerly referred M s = 980 ± 20 MeV and M P = 1300 ± 50 MeV HQ- F = 
90 ± 2 MeV [28, 50J and the standard reference ;^PT renormalization scale fiQ = 770 MeV. The 
short-distance constraints determine c m and d m in terms of the scalar and pseudo-scalar masses, 
producing 

L 8 = (0.83 ±0.05) • 10~ 3 , C 38 = (8.4 ± 1.0) • 10~ 6 . (112) 

Naively, if the uncertainty on the saturation scale is estimated by observing the variation with jj, 
in the range 0.5-1 GeV, one would expect the former values to be deviated from the actual ones 
at the order of AL 8 ~ 0.3 • lO" 3 , AC 38 ~ 5 • 10~ 6 . 

In order to go beyond the naive estimate of the subleading 1/Nc uncertainty, we consider 
now the one-loop contributions computed in previous sections. In a first approach, we consider 
just operators in the lagrangian with at most one resonance field [8]. At one-loop, in addition 
to the tree-level exchanges, one has the two-meson absorptive channels tttt, Vtt, Stt and Pit, 
determined by the scalar parameters c m and Cd, the pseudo-scalar coupling d m and the vector ones 
Gy and My. If we work in the WSR-renormalization-scheme for c m and d m , the short-distance 

constraints produce at NLO the same structure found from the large- Nc WSR, & m = ^ j^rzfip 
and = 4r ,>2 Ms » A /f 2 ■ The other three resonance parameters (cd, Gy,My) are fixed by means of 

^ 1X1 p — Ivl Q 

the logarithmic OPE constraints fl83|) . olq = = <m = 0, giving 

c d = 60±4MeV, G v = 93 ± 5MeV , My = 853 ± 28 MeV. (113) 

These numbers are found to be quite off the physical ones, q ~ 30 MeV, Gy ~ 60 MeV, My w 
770 MeV [H El [231 12H [2SI [271 12H1 EJJ. The LEC prediction for the standard comparison scale 
/io = 770 MeV then result, 

L 8 (/i ) = (2.28 ± 0.19) ■ 10" 3 , C 38 (/i ) = (26 ± 4) ■ HT 6 . (114) 
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In order to get these SU(3) %PT couplings, we employed in the U(3)-SU(3) matching Eq. (jlOOp the 
chiral singlet pseudoscalar mass m = 850 ± 50 MeV [50] . These estimates are still far from former 
values in the bibliography for /i = 770 MeV: L 8 = 0.9 • 10" 3 and C 38 = 10 • 10" 6 from 0(p G ) xPT 
and resonance estimates [5], later refined into L s = (0.61 ± 0.20) • 10~ 3 [6] and recently updated 
into L 8 = (0.37±0.17)-10" 3 [7]; L 8 = (0.6±0.4) • 10" 3 and C 38 = (2±6)-10" 6 from a previous NLO 
calculation in R^T [3T]; L 8 = (1.02±0.06) • 10~ 3 and C 38 = (3.3±0.6) ■ 10" 6 from Dyson-Schwinger 
equation analysis [52]; L 8 = (0.36 ± 0.05 ± 0.07) ■ 10~ 3 from Lattice simulations [53] . 

Although the calculation with just the Cr operators is able produce an appropriate description 
of the tttt channel (thanks to the Cd(Su il u tJ ') operator), its coupling q gets a extremely shifted 
value as this parameter has been used to accommodate the OPE at short distances. This does 
not represent by itself an important drawback in our analysis, where the goals are the LECs and 
R%T is devised as a convenient interpolator between high and low energies. However, the problem 
in our case is the erroneous description that one obtains for the Rtt channels with only the Lr 
operators [321 The Stt contribution to the spectral function behaves like a constant and the 
Vtt one grows with the energy, moreover, as My is also determined from the OPE matching, the 
position of the first two-meson threshold after the tttt one (i.e., the Vtt channel) is shifted from 
its physical place. 

8.2 Improving one Rtt channel: extending the lagrangian 

The straight forward procedure to ameliorate our one-loop amplitude is the inclusion of the re- 
quired operators for the proper description of the lowest absorptive cuts, this is, tttt and Vtt. The 
first one is ruled by the already included q operator but the latter demands the Af v term from 
Eq. (I102p . which now induces PVir interactions and allows to cure the infinitely growing behaviour 
of the Vtt contribution to the spectral function. 

Now we use the former inputs Ms, Mp, F, mo and the physical coupling a = 30 ± 10 MeV [H 
EBl EZl EEl EH] • The remaining parameters (Gy, My, Af v ) are extracted from the three logarithmic 
OPE constraints = — ot± = 0. Indeed, this system only has real solutions in the very 
corner of the parameter space, for low pseudo-scalar mass (Mp w 1.25 GeV) and high q and scalar 
mass (cd ~ 40 MeV, M$ ~ 1.00 GeV). This does not improve the value of the vector coupling and 
mass with respect to the former section, which become Gy ~ 120 MeV and My w 400 MeV. The 
LEC predictions result, 

L 8 (fi ) « 0.5 • KT 3 , C 38 M ~ -8 • 10" 6 , (115) 

where L 8 may look acceptable but the presence of such a low distorted Vn threshold is reflected 
in a value of C 38 which looks still a bit off. Nonetheless, these values are closer to those formerly 
obtained in the bibliography [5J (6J EH EH E3] . 

The problem is that the Vtt is not the only relevant channel that appears after the ittt one. The 
Stt channel opens up at an energy not far from the Vtt threshold. Thus, even if the Vtt channel 
can be now correctly described, the Stt contribution to the spectral function still shows a wrong 
constant behaviour [321 H6] . The Af p operator in (11021) is then crucial to cure that behaviour. 
Furthermore, this operator mends as well the similar bad short-distance behaviour found in the 
Pit cut contribution to the SS spectral function. 

Nonetheless, the presence of Af v in the lagrangian is still essential. If one repeats the NLO 
computation adding only the Af p operator (but not Af y ) the vector parameters become of the 
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Figure 14: Comparison of the LEC predictions in this work with previous results in the bibliography. 

order of G v ~ 20 MeV and My ~ 2 GeV. On the other hand, the LEC predictions L 8 ~ 1.3 • 1CT 3 
and C38 ~ 12 ■ 10~ 6 seem to improve with respect to the case with only Cr operators in the R%T 
lagrangian [8], with at most one resonance field. 

The inclusion of the \f A operator alone seems to move the results also in the right direction. 
Although it does not affect the previous channels, it opens the Air absorptive cut. Even if its effect 
at low energies is small, it helps to fulfill the OPE constraints. Taking now the extra needed input 
My = 770 ± 20 MeV together with the former ones, it is possible to extract the remaining ones 
(\f A ,Gy, Ma) through the three log OPE conditions. The value for the vector coupling turns 
out to be now more natural (Gy = 67 ± 18 MeV) but the ai(1230) mass falls down to very low 
values (Ma = 610 ± 50 MeV). The predictions for the chiral couplings show a clear improvement, 
L 8 = (0.7 ± 0.4) • 10~ 3 , C 38 = (4 ± 5) • 10~ 6 . 

8.3 Improving the Vir, Sir, An and Pit channels 

In order to have a proper description of all the Rti absorptive cuts, the Af A , \ f v and Af p operators 
from Eq. ( 11021) are now included in the R%T action. We take the same inputs as before, Ms = 
980 ± 20 MeV, M P = 1300 ± 50 MeV, F = 90 ± 2 MeV, m = 850 ± 50 MeV, c d = 30 ± 10 MeV, 
M A = 1230 ± 200 MeV, My = 770 ± 20 MeV and Gy = 60 ± 20 MeV. Both c d and Gy have 
been taken with a naive 33% error, as they appear only in the NLO part of the correlator. This 
will account for the possible NNLO variations in the one-loop correlator depending on whether 
it is evaluated with these physical couplings or their large-iVc values. The remaining unknown 
parameters (Af v , Af p , Af A ) are extracted from the three logarithmic OPE constraints, leading to 
our final LEC estimates, 

Lsifio) = (1.0 ±0.4) ■ HT 3 , C 38 (/io) = (8±5) ■ 10- 6 . (116) 

These numbers are compared to previous determinations in Fig. [141 Although there is still a clear 
dispersion between the various measurements, at the present error level we remain essentially 
compatible. Further efforts should be focused on the extraction of the scalar and pseudo-scalar 
pole masses in order to sizably reduce the uncertainties in the R/\T calculations. 

In general, the three logarithmic OPE constraints = — = produce complex 
solutions for the Af p , Af y , Af" 4 . In order to remain within the quantum field theory description, 
only the real values are kept. The regions with at least one real solution are shown in Fig. [15j 
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Figure 15: Allowed regions with one (light blue) or two (white) real solutions from the logarithmic OPE 
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constraints a = = a\ = 0. No real solution exists in the darker purple regions. In the upper left 
corner one can see the dashed rectangle provided by the ranges Cd = 30 ± 10 MeV, Gy = 60 ± 20 MeV. 

There, we plot the allowed ranges for a and Gy, with the other inputs taken at their central 
values. Indeed, there is no real solution for the central values q = 30 MeV and Gy = 60 MeV. 
On the contrary to other phenomenological analysis which seem to prefer a q coupling below 
30 MeV [23J E7J [28] , the log OPE constraints require slightly larger values, cd £ 30 MeV. However, 
in general for a around 30 MeV is always impossible to have real solutions for the value of the 
coupling Gy ~ 64 MeV obtained from V decays [8], [9j [28]. Actually, if one demanded the Tin 
scalar form-factor (and the corresponding tttt contribution to the SS spectral function) to vanish 
at high energies one would obtain Cd = F 2 /4c m ~ 42 MeV. However, in this work we do not 
perform a channel by channel analysis as in Ref. [31J. Indeed, in our field theory approach one 
could fix separately the short- distance behaviour of the tttt and all the Rtt channels through the 
X RR operators, but the latter also generate RR' absorptive cuts with the wrong properties at 
high momentum. The only option is the global adjustment of parameters considered in this work, 
where the lowest channels arrange the short-distance behaviour of the highest cuts at the price of 
slight modifications on their couplings. 

The allowed (c<j, Gy) region of Fig. [TS] actually changes if one varies the other inputs. Thus, 
we observed the whole range of the LECs allowed for the possible variations of the inputs and 
used this interval as our estimate of the central value and error. The maximum (minimum) value 
of the LECs was obtained at the largest (smallest) Cd and Gy. Likewise, the most extreme LEC 
values were obtain when Mp and Ma became smaller and M$ larger. These three parameters are 
responsible for most of the uncertainties. The impact of the My, F and mo errors in the global 
precision is negligible. 

The R/\T computation progressively approaches the physical value as one incorporates more 
and more physical information. This is quite non-trivial, as the introduction of a new chiral 
invariant operator leads to the opening of the new absorptive cuts in addition to those channels 
we are in principle interested in. For instance, the c m (Sx+) rules the decay into one scalar 
resonance and also contributes to the S'-meson exchange in the tttt channel. But at the same time 
it also induces the decay into Stt (though other operators like Af p are also relevant). Thus, the 
Crr' terms were used in our calculation to improve the description of the Rtt channels, which were 
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incompletely described by the linear lagrangian Cr [8]. The price to pay was that new RR' channels 
with two intermediate resonances showed up in our NLO computation of the correlator. Although 
the impact of these higher thresholds is suppressed at low energies if one chooses a convenient 
renormalization scheme [32, 46J, their impact in the high-energy matching and OPE constraints 
is a priori non-trivial. In this paper we find that, indeed, the most relevant information in order 
to extract the low energy chiral couplings seems to be provided by the lightest cuts. On the other 
hand, one realizes that the values of the couplings differ from those in the full large- Nc theory [40] 
and that the description of the heaviest absorptive channels may be very distorted [43J. Indeed, 
we obtain the resonance couplings Xf p = -0.22 ±0.08, Af y = 0.14 ±0.07 and \Xf A \ = 0.16 ±0.14. 
Even though these numbers have the right signs and order of the magnitude as the theoretical 

expectations Af p = = Cd 2 ~^ Cm 0.7, Af y = ^% 0.7 and Xf A = (in our analysis, for 

convention, we have took c m , d m and Gy as positive), their values are still far from being accurate 
determinations of these parameters. 

8.4 Impact of the RR' channels 

In this section we will make a digression on the importance of the RR' intermediate cuts that 
are opened after including the Crr> operators in the LO action. We will remove by hand the 
contributions with two-resonance cuts. Although this procedure is not well justified from the 
QFT point of view, we will perform this exercise in order make a rough comparison with the 
previous dispersive calculation of the octet SS — PP correlator [3T]. The RR' channels were 
neglected there, as their contribution in the dispersive integral was suppressed at low energies by 
inverse powers of {Mr + Mr*) 2 . 

Thus, we redid the calculation and removed by hand the diagrams with two-resonance cuts. 
This expression was then matched to the OPE at short distances, producing finally the low-energy 
constants, 

LsGuo) = (0.1 ± 0.7) ■ HT 3 , C 38 (/i ) = (-3 ± 9) • 10" 6 , (117) 

where we used the same inputs as in the previous subsection. The errors are now found to be 
larger and, though compatible with our final result (11161) . the elimination of the RR' cuts decreases 
slightly the range for the LEC determinations, approaching them to the lower values preferred by 
recent 0(p 6 ) analysis [7] and lattice simulations [53]. However, discarding these heavier channels 
from the one-loop computation in this way does not seem very sound from the theoretical point 
of view and it is shown here just as an exercise. 

9 Conclusions 

In this paper, we have performed the one loop QFT calculation of the two-point SS—PP correlator 
within R%T. We started with Ecker et al.'s lagrangian [5], containing only operators with at 
most one resonance field, and renormalized step by step all the relevant vertex-functions and 
propagators. Then we imposed OPE constraints on the full one-loop correlator, not on separate 
individual channels as it was performed in a previous NLO calculation [31]. Likewise, no short- 
distance constraint from other observables [39J was used in the present article. 
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After fixing part of our R%T couplings through these high-energy conditions, we expanded our 
result at low energies. Due to the chiral invariant structure of R%T, we were able to match the 
chiral logarithms and found predictions for the coupling constants L 8 (//) and C 38 (/^). The 
large discrepancy of these first numerical determinations with respect to the numbers found in the 
literature indicated that the simple Lagrangian Cr (with operators with at most one resonance 
field [S]) pointed out the need for a more complicated structure of the R^T action. The Cr terms 
could not fully describe the dynamics of all the two- meson intermediate channels: just the im 
channel description was adequately provided by the operators with at most one resonance field; 
all other channels (Vn, Sit. . . ) did not have the right short-distance behavior. Thus, beyond any 
numerical discrepancy in the LECs, the absence of operators with two an three resonance fields 
produces a severe theoretical issue at high energies [30] . 

In order to arrange the Rn cuts with one resonance and one Goldstone we add all the operators 
£rr' with two resonance fields relevant for the SS — PP correlator to the leading R%T lagrangian. 
These are the Af p , \± v and Xf A terms given in Eq. (11021) . The introduction of these operators 
produce a dramatic improvement. When only one of them is added to the action, the LEC 
predictions move in the right direction, i.e., towards the range of values found in previous studies. 
After considering all the three C rr> operators, we obtain the final values for /io = 770 MeV, 

L s (jM>) = (1-0 ± 0.4) • 10" 3 , C 38 M = (8 ± 5) • HT 6 , (118) 

in reasonable agreement with the values obtained through other approaches [5], [7J [311 [521 [53] . 
We want to remark, that this result is progressively approached as more and more complicated 
operators are added to the hadronic action. The terms of the lagrangian that rule the lightest 
channels result crucial and, thus, those determining heavier cuts not included in the analysis are 
expected to produce little influence. 

The essential difference with the previous dispersive calculation of the SS — PP correlator at 
NLO [31 J is the presence of RR' cuts in the present work. These intermediate channels automat- 
ically show up at the very moment we place the Crri operators in the R%T action. Although it 
is possible to demonstrate that the contribution from these heavy RR' cuts is suppressed at low 
energies [32, 46J, their impact in high-energy conditions such as the NLO Weinber sum-rules is 
pretty non-trivial. The difference between the present article and Ref. [31] could be taken as a 
crude estimate of the impact of neglecting those higher channels. 

In addition to the estimation of LECs, we also discussed some general issues about renormal- 
ization schemes within R%T. The use of the running MS masses Mr(h) was not very convenient 
as their meaning changed as one added new operators to the R%T action. Thus, they were re- 
expressed in terms of pole masses Mr. Likewise, we found that, with respect to the large-iVc 
WSR, the NLO Weinberg sum-rules (!86|) led to large uncertainties and variations for the values 
of c m and d m derived from them in the MS'-scheme. A more convenient subtraction scheme was 
found to minimize these uncertainties that stemmed from the high-energy matching whereas, on 
the other hand, it was found to leave the low energy prediction ([98]) unchanged (except for the 
improved accuracy in the resonance coupling determination from short- distance constraints). 
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A Running of the renormalized parameters with Cq + C>r 

When only operators with at most one resonance fields are considered in the RxT action [8], one 
finds before performing the meson field redefinition the running, 



<91n/i 2 512tt 2 F 2 



[16(4 - d 2 J — F 2 — lQc d c m + 4c 2 , + 2G 2 V ] , 



dM 2 s dM 2 P dc m dd r , 



ain/i 2 din /i 2 <91n/i 2 <91n/i 2 

8X S 3c 2 dX P 



din /i 2 16vr 2 F 4 ' din fx 2 ' 

dX "s = 3c d d\? 3 

dln/i 2 64tt 2 F 2 ' dlnfx 2 K ' 

After the renormalization one may then consider a convenient field redefinition that removes 
precisely the renormalized X$,p, and Af 8 . They (and their running) seem to disappear from 
the theory although their information is actually encoded in the renormalized effective couplings 
that remain in the action. Their running turns out to be then 

dL cS 3 

dhT? = ^^[ 16c -- 16 ^- F2 - 8c ^ + 2c ^ + 2 ^- 16c ^] ' 

dMf 2 3c 2 M| dMf 2 _ 

<91n/i 2 ~ 16vr 2 F 4 ' d\n^ 2 ' 

atf - ^ {iCdCm - F >), *e =0 . (120) 







<91n/i 2 64tt 2 F 4V " ain/i 

B On- shell scheme for c m and d m 

This would be a continuation of the pole-mass scheme. In addition to this, the renormalized 
on-shell couplings c m and d m are prescribed, respectively, by the real part of the residue of the 
correlator at the scalar and the pseudoscalar resonance poles [3"Tl I3"2"] . This was the scheme 
considered in the dispersive approach from Refs. [3TI [32] . The shift An with respect to the 
MS'-subtraction prescription is given up to NLO in 1/Nq by 

2c m Ac m = c 2 m -c 2 m = ^rRe$: s (M 2 s ) u - c 2 m ReZ r s ' (M 2 S ) U , 
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2 d m /\d m 



ct-d 



Re% P (M 2 P ) U - d 2 m ReE P ' (M 2 ^ . 



2\U 



m m 



2Bn 



In the case where only C p interactions are considered, one has 



2c m Ac 
2 d n i^d rn 



Ac d c m 3M 2 



F 2 128vr 2 



0. 



In- 



m 2 



C Feynman integrals 



The scalar integrals are 



A (M 2 ) 



dk° 



t(27i) d k 2 - M 2 + ie 
dk d 1 



B (p 2 ,M 2 ,M 2 ) , i ( 27r )d( jfe2 _ M 2 + ie )[(p_ A .)2_ M 2 + ie ] 



Using the formula in [30] we use the following expansions 
A (M 2 ) 



—M 2 I M 2 
16^ Aoo + ln ^ 



5 (P 2 ,0,0) 
5 (p 2 ,0,M 2 ) 
B (p 2 ,M 2 ,M 2 ) 
J(p 2 ,MlM 2 ) 



1 



Ano - 1 + In 



-p 



167T 2 | \ /i^ 

I I M 2 ( M' 2 \ ( ,r 

16^ A - + ln ^- 1+ ( 1 -^J ln ( 1 -M^ 



1 I M 2 

< + In — - 1 + a M In 



\<7 M - 1/ 



327T 2 



2 + 



M 2 -M 2 M 2 + M 2 



P 



M 2 a - M 2 



In 



M 2 



A 1/2 (j> 2 , M 2 , M 2 ) ^ ( [p 2 + \ 1/2 (p 2 , M 2 a , M 2 )] - {M 2 a - M x 



2^2 



V 



[p 2 -X 1 / 2 (p 2 } M 2 a} Ml)Y-(Ml-M ( 



2\2 

b . 



where a M = a/1 — 4M 2 /p 2 and X(x, y,z) = (x — y — z) 2 — Ayz. 



D Useful expansions 

Using expansions for x — > oo 



2x 



2x 2 
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21n(— x) ln(— x) 3 
1>{x) = -]n(- x ) + —^--±j r ±- — + 

;i-I)ln(l-x) = ln(-x)-^t^-- + -^ + 

x xx 2x z 
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